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DETECTION  OF  A  CHI-SQUARE  FLUCTUATING  TARGET 

IN  GAUSSIAN  NOISE 


INTRODUCTION 

Detection  of  a  target  in  an  echo-location  system  is  often  accomplished  through  the  use  of  a 
matched- filter  receiver.  Shown  in  Fig.  1,  it  is  composed  of  a  multiplier  to  form  the  product  of 
the  return  (echo)  x{t)  and  the  processing  signaJ  g{t),  which  is  followed  by  an  integrator,  and  an 
envelope  detector.  The  output  r  is  then  compared  to  a  threshold  rj  to  make  one  of  two  decisions; 

Hq  (no  target  present):  x(t)  =  n(t),  (1) 

where  n{t)  is  white  Gaussian  noise,  or 

Hi  (target  present): i(<)  =  pe^f{t  —  t)  +  n(t),  (2) 

where  f{t)  is  the  transmit  signzd,  r  is  the  r£inge-delay  of  the  target,  and  p  >  0  and  G  are  the 
range  and  phase  of  the  reflection  coefficient  associated  with  the  point  target.  The  target  strength 
is  w  =  This  represents  a  fairly  simple  model  of  the  echo- location  problem;  a  more  complicated 
model  of  the  return  would  include  the  effects  of  target  Doppler  or  range-spreading.  The  output  of 
the  matched-filter  is  given  by 

rir)  =  |y  ^  dt  ,  (3) 

where  f  is  the  hypothesized  range-delay  of  the  t^get.  Note  that  g{t)  is  taken  to  be  a  delayed  and 
normalized  replica  of  the  transmit  signal  /(t).  For  our  purpose,  we  eissume  that  we  have  correctly 
hypothesized  the  range  delay  of  the  target,  that  is,  f  =  r.  This  implies  that  r  =  r(T)  =  [pe^  +  nj, 
where  n  is  the  matched-filter  response  to  the  noise. 


x(t) 


g*(t) 

Fig.  1  -  The  matched-filter  receiver  structure 

The  task  is  to  find  the  probability  distribution  and  cumulative  probability  function  p(r)T+c 
given  that  the  noise  is  a  white  Gaussian  noise  process,  and  that  p^  is  a  chi-square  random  variable  ^ 
This  will  allow  us  to  quantify  the  performance  of  the  matched-filter  for  target  detection. 

Manuscript  approved  November  16,  1993. 

*T,  C,  and  T  +  C  are  used  to  denote  the  density  function  and  cumulative  probability  function  of  the  amplitude 
of  the  matched-filter  response  to  the  target,  clutter,  and  target  plus  clutter,  respectively. 
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It  should  be  noted  that  an  expression  for  the  density  function  for  the  matched-filter  output  r 
under  the  assumptions  given  here  has  been  reported  by  others  [1,  2],  A  short  derivation  of  the  same 
result  is  given  here  with  the  intent  of  producing  an  available  reference.  Furthermore,  we  derive 
expressions  for  the  cumulative  probability  function  and  moments  of  the  matched- filter  output.  To 
the  best  of  the  author’s  knowledge,  these  results  have  not  been  reported. 

DERIVATION  OF  THE  DENSITY  FUNCTION 

Consider  the  case  of  using  a  matched-filter  for  detecting  a  nonfluctuating  target  (p  fixed)  in 
white  Gaussian  noise.  The  probability  density  function  for  r  is  the  well  known  Rician  density 
function  [3-5]  given  by 

p{Ap)t+c  -  '  (4) 

where  Iq  denotes  the  modified  Bessel  function  of  order  zero  and  is  the  expected  value  of  when 
only  noise  is  present  in  the  return^.  This  density  function  can  be  derived  assuming  that  the  target 
phase  0  is  either  fixed,  or  is  a  uniform  random  variable  that  assumes  its  values  in  the  interval 
[0, 27r).  Thus,  explicit  use  of  the  target  phase  in  the  derivation  that  follows  is  unnecessary,  and  will 
be  eliminated. 

The  target  strength  w  is  often  modelled  as  a  chi-square  random  variable,  for  which  the  distri¬ 
bution  is 

(5) 


p{w)t  = 


nN 


where  N  is  any  positive  real  number.  However,  it  is  more  convenient  to  use  the  density  function  of 
the  tcirget  reflection  amplitude  tt-  Thus,  by  the  transformation  of  random  variable  rj  =  y/w  [6], 
the  density  function  for  rj  is 


firr)  = 


2N^ 


riN)a 


2n‘''t 


(6) 


for  which  it  can  be  shown  that 


E{p}  = 


T{N  + 1/2) 
nN) 


E{p^}  =  <7. 


(7) 

(8) 


Note  that  Eq.  (8)  shows  the  (expected)  target  strength.  This  distribution  defines  what  is  known 
as  the  Nakagami  target  fluctuation  model  [2]. 

Figure  2  shows  the  curves  of  the  density  function  in  Eq.  (6)  for  several  values  of  N .  Note  that 
as  N  increases,  the  variamce  tends  to  zero,  and  the  density  function  approaches  a  delta  function. 
Thus,  as  N  increases,  the  target  amplitude  fluctuations  decrease,  and  in  the  limit  as  N  oo  the 
target  is  nonfluctuating.  Moreover,  E{p}  — >  y/a  as  N  oo. 

*Since  r  >  0,  p  >  0,  and  w  >  0,  it  should  be  understood  that  the  density  functions  and  cumulative  probability 
functions  given  in  this  report  are  only  defined  for  positive  semi-definite  values  of  their  arguments. 
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Fig.  2  -  Probability  distributions  of  the  target  reflection  coefRcient  amplitude 

For  certain  values  of  N,  the  density  function  in  Eq.  (6)  reduces  to  well-known  target  fluctuation 
models;  N  =  I  corresponds  to  a  Rayleigh  fluctuating  target,  and  N  =  2  corresponds  to  a  one- 
dominant-plus-Rayleigh  fluctuating  target. 

Given  p{r\p)j'+c  and  p(p)t,  it  follows  that 

/■OO 

Pir)T+C=  /  p{r\p)T+CPip)Tdp.  (9) 

JQ 

Thus,  substituting  Eq.  (4)  and  Eq.  (6)  in  Eq.  (9),  using  the  formula  ([7],  page  716,  Eq.  6.631.1), 


f 

JO 


(2  ■^2’^’ 


4a 


(10) 


where  Jo  denotes  the  Bessel  function  of  order  zero,  and  iFi  denotes  the  confluent  hypergeometric 
function,  and  realizing  that  Jo(jx)  =  /o(i),  one  finds  that 


(11) 


Two  other  forms  of  the  density  function  in  Eq.  (11)  can  be  found.  If  we  define  the  signal-to-noise 
ratio  (SNR)  as 


B{r^) 


1^’ 

then  the  density  function  in  Eq.  (11)  can  be  rewritten  as 

2re'’'*/*'" 


p(r)T+C  = 


62 


- i^iFi  iV;  1;  — ^ 

(i+f)"  V 


/62 


N 

77/ 


(12) 


(13) 


Furthermore,  if  we  define  y  to  be  the  matched-filter  output  normalized  by  the  rms  value  of  the 
(un-normalized)  matched-filter  response  to  noise,  i.e.. 
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y  = 


r 

v 


then,  by  the  linear  tramsformation  of  a  random  variable  [6], 


p(y)T+c  = 


2ye-y^ 


N 


iFi 


(14) 


(15) 


Calculation  of  the  density  function  is  easily  accomplished  using  Mathematical  a  software  pack¬ 
age  for  performing  numerical  and  symbolic  mathematics  [8],  since  the  hypergeometric  functions  are 
available  as  function  calls.  A  user-defined  Afathemattca  function  for  the  density  function  in  Eq.  (15) 
is 


DensityCn_,R_,y_]  :■  2  y  Exp[-y“2] 

HypergeometriclFl [n, 1 , (y‘2)/(l+n/R)] / (1+R/n) 'n ; 

Figures  3  and  4  show  plots  of  the  probability  density  functions  for  R  equal  to  1  and  10,  zind  N 
equal  to  0.5,  1,  3,  and  7.  Note  that  for  i?  =  1  the  density  functions  do  not  differ  significantly.  This 
is  reasonable  since  it  can  be  shown  that  in  the  limit  as  i?  0,  the  density  function  in  Eq.  (15) 
becomes  equal  to  2ye~*'  . 


DERIVATION  OF  THE  CUMULATIVE  PROBABILITY  FUNCTION 


We  note  that  the  confluent  hypergeometric  function  has  the  series  representation 


^0 


(16) 


where  (a),!  =  a(a  -b  1) . . .  (a  -f  n  —  1),  and  (a)o  =  1,  ^lnd  is  known  as  Pochhammer’s  symbol. 
Substituting  Eq.  (16)  into  Eq.  (15),  integrating  from  zero  to  y,  using  the  definition  of  the  incomplete 
gamma  function  defined  as 
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(17) 


and  realizing  that  (l)n  =  «!,  one  finds  that  the  cumulative  probability  function  is  given  by 


(18) 


As  with  the  probability  density  function,  calculation  of  the  series  in  Eq.  (18)  is  easy  since  the 
incomplete  gamma  function  is  available  as  a  function  call  in  Mathematica.  Of  course  a  truncated 
form  of  the  series  must  be  used,  and  it  was  found  that  the  terms  in  the  series  often  increased  before 
they  decayed.  Thus,  practical  experience  revealed  that  the  first  10  terms  should  always  be  used 
when  summing  the  series.  Once  beyond  the  tenth  term,  the  summation  can  be  terminated  if  the 
curren',  term  differs  from  the  previous  term  by  less  than  0.000001.  A  user-defined  Mathematica 
function  for  the  calculation  of  the  cumulative  probability  function  in  Eq.  (18)  using  the  method 
described  here  is 


CumProb Cn_ , R_ , y_]  :=  Module [{error,  p,  plast,  k,  temp}, 

error  =l;p=0;k*-l; 

While[(error  >  0.000001)  ||  (k  <  10) , 
plast  =  p; 
k  =  k+1; 

temp  =  Pochhammer[n,k]  Gamma[k+1 ,0,y*2] 
/  (Factorial [k] “2  (l+n/R)*k): 
p  =  p  +  N[temp/(l+R/n) “n] ; 
error  =  Abs [p-plast] ] ; 
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Figures  5  and  6  show  the  cumulative  probability  functions  for  R  equal  to  1  2ind  10,  and  N 
equal  to  0.5,  1,  3,  and  7.  An  expanded  horizontal  axis  was  used  in  Fig.  5  since  the  curves  are  nccU’Iy 
indistinguishable . 


Fig.  5  -  Cumulative  probability  function  for  Jl  =  1,  and  for  N  equal  to  0.5,  1,  3,  and  7 


Fig.  6  -  Cumulative  probability  functions  for  R  =  10,  and  for  N  equal  to  0.5,  1,  3,  and  7 

The  convergence  of  the  series  in  Eq.  (18)  can  be  proved  by  using  the  binomial  theorem.  First, 
we  note  that  from  Eq.  (17),  and  the  definition  of  the  gamma  function,  it  follows  that 

dx  =  r(n  +  1)  =  n!  . 


7(n+  l,y2)  <  f 

Jo 


(19) 
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This  inequality  implies  that 


PiyW+c  < 


oo 


E 


w, 

n! 


(20) 


It  can  be  shown  that  (N)n/n\  is  equal  to  the  binomial  coefficient.  Furthermore, 
1/(1  +  N/R)  <  1  since  R  >  0.  Thus,  by  the  Binomial  Theorem,  the  right  side  of  Eq.  (20) 
converges,  yielding 


■P(y)r+c  < 


1 


(i  +  f)  0  -  1+:^) 


N 


=  1 


(21) 


This  bound  implies  that  the  series  in  Eq.  (18)  converges. 


MOMENTS  ABOUT  ZERO 


The  m-th  moment  about  zero  is  defined  as 


=  /  y"'p(y)T+c«^y- 

Jo 

Substituting  Eq.  (15)  into  Eq.  (22),  and  using  the  formula  ([7],  page  860,  Eq.  7.621.4) 
jf  e"**x^~SFi(a;c;aa;)da:=  ^^2F\  b;  c;  , 
where  2F1  is  the  hypergeometric  function,  one  finds  that 


+  1;  1; 


1  +  ^ 

^  ^  R  . 


However,  the  hypergeometric  function  obeys  the  property 

7— 2F1  (a,  6;  c;  — ^  )  =  2^1(0,  c  -  b;  c;  -i), 

(l  +  i)“  1  + 

and  can  be  used  to  derive  ^ln  alternate  form  of  the  m-th  moment  about  zero  given  by 


(22) 


(23) 


(24) 


(25) 


(26) 


The  hypergeometric  function  has  the  series  representation 


OO 


2Fi(a,  b;  c;x)  =  ^ 

n=0 


(a)n(b)n  X" 

(c)„  nr 


(27) 


which  is  convergent  for  |a:|  <  1.  However,  if  the  series  terminates,  it  is  convergent  for  all  x.  One  can 
exploit  this  property  by  noting  that  for  m  =  2fc  (fc  =  0, 1,2, . . .)  ,  the  series  for  the  hypergeometric 
function  in  Eq.  (26)  terminates  at  the  fc-th  term,  thus 
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=  fc! 


k 

l  +  ]^fc(fc-l)...(jk-n  +  l) 

n=l 


{N)n 

(n!)2 


(28) 


In  p2irticul£U‘, 


£;{y2}  =  i  +  i?. 


(29) 

(30) 


Of  course,  Eq.  (28)  can  be  used  to  prove  the  trivial  result  E{y®}  =  jB{1}  =  1. 

A  user-defined  Mathematica  function  that  calculates  the  moments  about  zero  is 
Moment  Cn_ , R_ , m_]  : ■ 

Gamma [m/2+1]  Hypergeometric2Fl [n, m/2+1 , 1 , 1/ (1+n/R)] /(1+R/n) ‘n; 


CALCULATION  OF  DETECTION  PROBABILITIES 


If  no  target  is  present,  then  i2  =  0,  and,  as  indicated  in  Section  2,  Eq.  (15)  reduces  to 

P(y)c  =  2ye"S'*.  (31) 

If  the  detection  threshold  is  rj,  then  the  probability  of  false  alarm  is 

P/a  =  f  2ye~*'*  dy  =  e"''*.  (32) 

Jri 

Solving  for  the  threshold  in  Eq.  (32)  yields 

V  =  y/-  ln(P/a).  (33) 

The  associated  probability  of  detection  follows  from  Eq.  (18): 

/•oo 

Pd=  f  p{y)T+cdy 

Jj] 

=  1-/  p{y)T+cdy 
Jo 

^  {N)nj(n  +  l,Ti^) 

“  +  (!  +  «)”  • 


Figmes  7  through  9  show  the  probability  of  detection  for  N  equal  to  0.5,  1,  3,  and  7,  and  a 
probability  of  false  alarm  of  0.01,  0.001,  and  0.00001. 
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CONCLUSIONS 

We  have  presented  a  concise  derivation  of  the  probability  density  function,  cumulative  proba¬ 
bility  function,  and  moments  about  zero  of  the  normalized  matched-filter  output  when  detecting  a 
target  buried  in  white  Gaussiam  noise,  amd  whose  tairget  strength  follows  a  chi-squaue  distribution. 
It  was  shown  that  the  cumulative  probability  function  has  a  series  representation  that  lends  itself 
te  calculation  using  a  mathematical  software  package  such  as  Mathematica. 


Fig.  7  -  Probability  of  detection  as  a  function  of  SNR  for  Pf,  =  0.01 


Fig.  8  -  Probability  of  detection  as  a  function  of  SNR  for  P/„  =  0.001 
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Fig.  9  -  Probability  of  detection  as  a  function  of  SNR  for  P/,  =  0.00001 
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